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Turbulent viscosity in cold accretion disks is likely to be hydrodynamic in origin. We investigate the growth 
of hydrodynamic perturbations in a small region of a disk, which we model as a linear shear flow with Coriolis 
force, between two parallel walls. Although there are no exponentially growing eigenmodes in this system, 
because of the non-normal nature of the modes, it is possible to have a large transient growth in the energy of 
certain perturbations. For a constant angular momentum disk, the energy grows by more than a factor of 1000 
for a Reynolds number of only 1000, and so turbulence is easily excited. For a Keplerian disk, the growth is 
more modest, and energy growth by a factor of 1000 requires a Reynolds number of nearly a million. Accretion 
disks have even larger Reynolds numbers than this. Therefore, transient growth of perturbations could seed 
turbulence in such disks. 



1. INTRODUCTION 

The origin of hydrodynamic turbulence is still not 
completely understood. Many efforts have been de- 
voted to solve this problem, beginning with the work 
of Kelvin, Rayleigh, and Reynolds at the end of the 
nineteenth century, to more recent work. In most 
cases, a significant mismatch has been found between 
the results based on linear theory and experimen- 
tal (observational) data. In the laboratory, plane 
Poiseuille and plane Couette flows become turbulent 
at Reynolds numbers, R > R c c xp ~ 1000, 350, re- 
spectively, while theoretically, plane Couette flow is 
linearly stable at all values of R and plane Poiseuille 
is linearly unstable only for R > Rf 10 = 5772. What 
is the reason for this mismatch? It would appear that 
for R values in between R = R e c xp and R — R^, 
some new effect allows the system to make a subcriti- 
cal transition to turbulence in the absence of a linear 
instability. 

In the field of astrophysics, hydrodynamics and tur- 
bulence find extensive applications. Accretion flows in 
binary stars, young stellar objects and active galactic 
nuclei must be turbulent in order for the gas to ac- 
crete. However the actual origin of the turbulence is 
unclear and still under debate. Inward mass accre- 
tion in a disk occurs as a result of the transfer of an- 
gular momentum outward by means of some viscous 
torque in the system. More than three decades ago, 
Shakura & Sunyaev and Lynden-Bell & Pringle 0] 
proposed that turbulent viscosity provides the neces- 
sary viscous torque to cause inward mass transport. 
However, the physical origin of the turbulence in an 
accretion disk was unclear until the work of Balbus 
& Hawley who showed that magnetized disks have 
a Magneto-Rotational-Instability (MRI) d H in the 
presence of a weak magnetic field. This instability 
provides a natural means to generate turbulence and 
transport angular momentum outward. Later, Haw- 
ley, Gammie & Balbus 6] showed that the turbulence 
dies out if the Lorentz forces are turned off in a magne- 
tohydrodynamically turbulent Keplerian disk. Subse- 



quently they again showed [7j that the magnetic field 
dies out when the tidal and Coriolis forces are turned 
off, keeping the Lorentz forces. In two subsequent 
papers [8|, Q , it was shown through numerical simu- 
lations that, whereas pure hydrodynamic turbulence 
is easily triggered in plane Couette flow and in a con- 
stant angular momentum disk, turbulence does not 
develop in an unmagnetized Keplerian disk even in 
the presence of large initial perturbations. The au- 
thors argued on this basis that hydrodynamic turbu- 
lence cannot contribute to viscosity in accretion disks. 

Despite the above results, it is probably prema- 
ture to rule out hydrodynamic turbulence in rotation- 
ally supported disks. Many astrophysical systems are 
known, e.g., star- forming disks, cataclysmic variables 
in quiescence, outer regions of AGN disks, etc., in 
which the gas is cold and neutral so that there is negli- 
gible coupling between the magnetic field and the gas 
(see e.g., [l(j[ll|,[l2|). ^ n these systems, the Magneto- 
Rotational-Instability (MRI) cannot play a role. How 
do these systems sustain mass transfer in the absence 
of the MRI? What drives their turbulence? It is pos- 
sible that the turbulence is purely hydrodynamic in 
origin. Even in the absence of exponentially growing 
modes, plane Couette flow is known to exhibit a large 
transient growth in energy for certain initial condi- 
tions. Is this transient growth responsible for the sub- 
critical turbulence that occurs in this flow, and could 
a similar phenomenon be responsible for generating 
hydrodynamic turbulence in cold accretion disks? Re- 
cent laboratory experiments on rotating Couette flow 
in the narrow gap limit with linearly stable rotational 
angular velocity profiles (similar to Keplerian disks) 
seem to indicate that turbulence does manage to de- 
velop in such flows . Longaretti Q points out that 
the absence of turbulence in the simulations [1,0 may 
be because of their small effective Reynolds number. 

Transient growth occurs in a linear system because 
of the non-normal nature of the associated operator 
pH Hil lrj| . How, and if, such growth might fi- 
nally cause turbulence is not well understood. One 
possibility is that the growth causes fluctuations to 
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acquire large amplitudes, and non-linear interactions 
then push the system into self-sustained chaotic tur- 
bulence. 

We present here some results on shearing rotating 
flows following the methods that have been developed 
to study the evolution of non-normal modes in plane 
Couette flow. We employ a local approximation, and 
focus mainly on an eigenvalue analysis of perturba- 
tions in Eulerian coordinates (for details, see [l^, and 
also |20] for a similar analysis in Lagrangian coor- 
dinates). In the next section, we briefly describe the 
basic model and write down the set of equations to 
be solved. Then in §3, we present a few important 
results. We conclude with a discussion in 84. 



C(x)e J<Tt e l ( kyV+kzZ \ where a = an + iaj is an eigen- 
value of the problem. After substituting these per- 
turbations in (JJJ and (J2J), we solve the combined Orr- 
Sommerfeld and Squire equations to obtain the eigen- 
system {a, u, £}. The detailed set of equations and 
its reduction to an eigenvalue equation form and the 
solution procedure are described in . 

As mentioned above, even in the absence of any ex- 
ponentially growing mode in the system, plane Cou- 
ette flow can exhibit a large transient growth in the 
energy of certain perturbations 0, • This growth 
occurs in the absence of non-linear effects and is be- 
lieved to facilitate the transition from laminar to tur- 
bulent flow. The maximum growth in the perturbed 
energy can be expressed as 



2. MODEL 

We consider a small patch of an accretion disk 
whose local geometry is described in terms of Carte- 
sian coordinates, x = [r — r ), y = r (cf) — (p ). We 
assume that the flow is incompressible and extends 
from x = — 1 to +1, between two rigid walls with no- 
slip boundary conditions. The y and z directions are 
unbounded with periodic boundary conditions. This 
flow behaves like rotating Couette flow in the narrow 
gap limit, when the unperturbed velocity corresponds 
to a linear shear, U = (0, —x, 0), and the Coriolis force 
is described by the angular frequency uj — (0, 0, 1/q): 
fl = rio(ro/r) q . Here all the coordinates and quanti- 
ties are expressed in terms of dimensionless variables. 
Note, q — 3/2 corresponds to a Keplerian disk, q = 1 
and 2 to a flat rotation curve and a constant angular 
momentum disk, respectively, and q = corresponds 
to rigid body rotation. A schematic diagram of the 
flow in the local coordinates is shown in Fig. 

Writing the perturbations as U x (x,y,z,t) — > 
u(x,y,z,t), U y (x,y,z,t) -> U y {x) + v(x,y,z,t), 
U z {x,y,z,t) -> w(x,y,z,t), P/p -> P/p+p(x,y,z,t), 
and combining the linearized Navier-Stokes and con- 
tinuity equations, we obtain 
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where £ = dw/dy — dv/dz. Equations Q and J2J 
are the standard Orr-Sommerfeld and Squire equa- 
tions, respectively, except that they now have addi- 
tional terms due to the presence of rotation in the 
system. We have to solve equations and @ with 
the boundary conditions: u — du/dx = £ = 0, at 
x = ±1 (no slip). 

Due to translational invariance in the y and 
z directions, the perturbations can be chosen as 
u(x,y,z,t) = u(x)e l,Jt e l{ - k y v+kz7 ^ and ((x,y,z,t) = 



where 



1 



optimum [ — 



,9.+ , du^ du ;+ j 
k 2 v)u+— 7r +C t C 
dx dx 



(3) 



dx (4) 



and Eq is the initial energy of the perturbation at time 
t = 0. By "optimum" we mean that we consider all 
possible initial perturbations and choose that function 
that maximizes the growth of energy at time t. To 
understand the technical details of how to compute 
the growth, see [l8|. 



3. RESULTS 

The eigenvalue analysis shows that, for plane Cou- 
ette flow and for a rotating flow with q < 2, there are 
no exponentially growing modes in the system, i.e., 
for no choice of the parameters is there an eigenvalue 
with 01 > 0. However, for plane Couette flow, previ- 
ous workers (e.g. [TEl IT^|') have shown that there is 
a large transient growth even for R = 1000. The in- 
teresting point is that the maximum growth in energy 
for plane Couette flow and a constant angular momen- 
tum disk (q — 2) are very similar. Figure [21 shows the 
contours of constant maximum growth (G max ) and the 
time at which the growth is maximized (i max ) for both 
plane Couette and q = 2. Whereas for q — 2, G max 
occurs exactly on the k z axis, for plane Couette flow 
it is slightly off the axis. This is the only difference 
between the two cases. 
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Figure 1: Unperturbed flow in the local shearing box. The arrows represent the velocity field. 
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Figure 2: Contours of (a) G max for plane Couette flow, where dashed contours correspond to 10,30, ...90; solid to 
150,200, ...600; dotted to 700,800, ...1100; dot-dashed to 1150, 1160, ...1200; (b) t max for plane Couette flow, where 
dot-dashed contours correspond to 20, 25, ...45; solid to 50, 55, ...100; dotted to 110, 120, ...180; (c) G max for q = 2 flow, 
where contour labels are the same as in (a); (d) £ ma x for q — 2 flow, where contour labels are the same as in (b). All 
calculations are for R = 1000. 



Table 1 

Maximum Growth for Various q and R — 1000 



Q ky k z 


Cm ax 


^max 


1.5 1.21 


13.04 


8.8 


1.7 1.06 0.74 


13.4 


9 


1.9 0.54 1.44 


23.75 


12.4 


1.99 0.04 1.8 


122.8 


27.3 


2 1.66 


1165.9 138.5 



Figure |3 shows how G max and i max vary as a func- 
tion of Reynolds number, R, for plane Couette flow 
and q = 2 flow. We see that G max scales as R 2 and 
t max as R in both cases, again highlighting the sim- 
ilarity of the two flows. The physical reason is that 
a constant angular momentum disk has zero epicyclic 
frequency. Thus, the effect of rotation is essentially 
cancelled out, making the basic structure of the sys- 
tem very similar to that of plane Couette flow. 

For a given R, when q slightly deviates from 2, 
growth immediately reduces dramatically due to the 
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Figure 3: Solid and dotted lines correspond to plane Couette and a q = 2 disk respectively. 



introduction of a small epicyclic frequency (= k = 
\/2(2 — q)fi) . This indicates the dominant effect that 
k has on the fluid dynamics. Table 1 and Fig. 0Jl 
show how G max and i max change as q decreases below 
2. It is interesting to note that, for a q = 2 flow the 
growth is maximum for k y = 0, while for a Keplerian 
disk (q = 1.5) the growth is maximum for k z = 0. 
Therefore, the location of G max moves in the k y — k z 
plane systematically from one axis to the other, as q is 
varied. Thus, for a constant angular momentum disk, 
we need to include vertical structure in the perturba- 
tions to understand the energy growth, whereas for a 
Keplerian disk a 2-dimensional analysis is sufficient. 

Even though the growth decreases dramatically as 
q falls below 2, nevertheless, at a given q, if the 
Reynolds number R increases to a large enough value, 
the growth can become significant. In Fig. we 
show how the peak value of growth increases with in- 
creasing R for Keplerian disks. It is found that G max 
scales as i? 2 / 3 and i max as i? 1 / 3 for large R. Therefore, 
while the presence of a finite epicyclic frequency has 
a strong dynamical effect on growth, it does not rule 
out a large growth. One just needs much larger values 
of R. Cold astrophysical accretion disks can have R as 
high as 10 10 or more, so large growth should in prin- 
ciple be easily achieved in such disks. Also at smaller 
k y , G max scales as k y , while at larger k y , G max de- 
creases as k y 4 ^ 3 . At large t, i max scales as k y 2 ^ 3 . For 
a detailed derivation of these scaling relations, see [Isj . 

For a Keplerian flow the maximum growth occurs 
for k y ~ 1.2, k z = 0. Figure shows the development 
with time of the perturbed velocity component u(x, y) 
corresponding to R — 1000. We show snapshots corre- 
sponding to four times: t = 0, t max /2, t roax , 3i max /2. 
The perturbations are seen to resemble plane waves 
that are frozen in the shearing flow. The initial per- 



turbation at t — is a leading wave with negative 
x-wavevector k x and with \k x \ k y . With time, the 
wavefronts are straightened out by the shear, until at 
t = frmaxj the wavefronts are almost radial and k x ~ 0. 
Then, at yet later times, the wave becomes trailing 
and the energy also decreases. The perturbations are 
very similar to the growing perturbation described by 
other authors [2ll I22. . 



4. DISCUSSION 

We find that significant transient growth of pertur- 
bations is possible in a shear flow with Coriolis force 
between walls (see also 24j). This system is an ideal- 
ized local analog of an accretion disk. Although the 
system does not have any unstable eigenmodes, never- 
theless, because of the non-normal nature of the eigen- 
modes a significant level of transient energy growth is 
possible for appropriate choice of initial conditions. 
If the maximum growth exceeds the threshold for in- 
ducing turbulence, it is plausible that this mechanism 
could drive the system to a turbulent state. Presum- 
ably, once the system becomes turbulent it can re- 
main turbulent as a result of nonlinear interactions 
and feedback among the perturbations. 

In this mechanism of turbulence, the maximum en- 
ergy growth and the time needed for this growth are 
probably the main factors that control the transition 
to hydrodynamic turbulence. As mentioned in §1, for 
plane Couette flow R e c xp ~ 350 and according to our 
analysis, for R = 350, G max = 145, and the corre- 
sponding t roax = 42.3. Since a q = 2 disk is very 
similar to plane Couette flow, the critical Reynolds 
number for turbulence in this case is also likely to be 
R c - 350. For this R, G max = 143.5 and i max = 48.3. 
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Figure 4: Growth curves for (a) q — 1.5 (solid), 1.7 (dotted), 1.9 (dashed), 1.99 (dot-dashed); other parameters are 
given in Table 1. (b) Growth curves for q = 1.5 when {R = 4000, k y = 1.2} (solid), {R = 1000, k y = 1.21} (dotted), 
{R = 500, k y = 1.29} (dashed); k z = 0. 
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Figure 5: Development of the perturbed ^-component of the velocity u(x, y) as a function of time for the optimal 
perturbation with the maximum growth of energy in a Keplerian flow with R = 1000. The perturbation has k y = 1.21, 
k z — 0, t miU ; = 8.8. Panels show the perturbation at (a) t = 0, (b) t = i m ax/2 = 4.4, (c) t = t max = 8.8, (d) 
t — 3i max /2 = 13.2. Solid and dotted contours correspond to positive and negative values of u respectively. 



Based on the above idea, we make the plausible 
assumption that the threshold energy growth factor 
needed for transition to turbulence in a shear flow 
is E c ~ 145. Applying this prescription to the opti- 
mal two-dimensional perturbations of a Keplerian disk 
(§§2,3), we estimate the critical Reynolds number for 
a Keplerian flow to be R c ~ 3.4 x 10 4 , i.e., a factor 



of 100 greater than in the case of plane Couette flow. 
The corresponding i max = 28.3, which is comparable 
to that in plane Couette flow, and is not too large 
compared to the accretion time-scale of a geometri- 
cally thin disk. 

Instead of taking R c ~ 350, we might wish to be 
conservative and take R c ~ 1000 for plane Couette 
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flow and q = 2 flow. In this case, G max ~ 1200 and 
imax ~ 120 — 140. Applying the requirement of E c ~ 
1200 to a Keplerian flow, we find R c ~ 10 6 and i max ~ 
100. Now the critical Reynolds number is a factor of 
1000 greater than in the case of plane Couette flow. 

An interesting result is that epicyclic motions in 
a differentially-rotating disk kill the growth dramat- 
ically and as a result the critical Reynolds number 
R c becomes higher. This also changes the optimum 
wavevector {k y ,k z } of the perturbations needed to 
produce energy growth. For a constant angular mo- 
mentum disk (q = 2) and plane Couette flow, it is 
seen that growth is maximized for k y ~ 0. Even for a 
very small shift in the value of q below 2, the location 
of maximum growth moves significantly in the k y — k z 
plane away from the k z axis. With decreasing q, the 
epicyclic motions of the disk increase, and correspond- 
ingly the optimum value of k y for growth increases 
while the optimum k z decreases. For a Keplerian disk, 
the growth is maximum for k z = (on the k y axis). 
To the best of our knowledge, this change in the loca- 
tion of the maximum growth in the k y — k z plane has 
not been commented upon prior to this work. 

In earlier numerical simulations 0, the above 
reduction of growth due to the effect of epicyclic mo- 
tion in the disk was already noticed. However, those 
authors then proceeded to rule out the possibility of 
hydrodynamic turbulence in Keplerian disks. We do 
not agree with this conclusion. As we have shown, 
Keplerian disks can support large transient energy 
growth, only they need much larger Reynolds numbers 
to achieve the same energy growth as plane Couette 
flow or a q — 2 flow. The numerical simulations prob- 
ably had effective Reynolds numbers < 10 4 (because 
of numerical viscosity) which is below our most op- 
timistic estimate of the critical Reynolds number for 
a Keplerian flow. Thus, we suspect the simulations 
simply did not have sufficient numerical resolution to 
permit the turbulence. The same point has been made 
by other authors [TiE^. 

We conclude with an important caveat. While the 
demonstration of large energy growth is an important 
step, it does not prove that Keplerian disks will neces- 
sarily become hydrodynamically turbulent. Umurhan 
& Regev |23| have shown via two-dimensional sim- 
ulations that chaotic motions can persist for a time 
much longer than the time scale t ma x needed for linear 
growth. However, they also note that their perturba- 
tions must ultimately decline to zero in the presence 
of viscosity. To overcome this limitation, it is neces- 
sary to invoke three-dimensional effects. Secondary 
instabilities of various kinds, such as the elliptical in- 
stability (e.g. [23), are widely discussed as a possi- 
ble route to self-sustained turbulence in linearly per- 



turbed shear flows. It remains to be seen if these in- 
stabilities are present in perturbed flows such as those 
shown in Figure 
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